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Abstract 

Cilia and flagella are hair-like extensions of eukaryotic cells which generate oscillatory beat 
patterns that can propel micro-organisms and create fluid flows near cellular surfaces. The evo- 
lutionary highly conserved core of cilia and flagella consists of a cylindrical arrangement of nine 
microtubule doublets, called the axoneme. The axoneme is an actively bending structure whose 
motility results from the action of dynein motor proteins cross-linking microtubule doublets and 
generating stresses that induce bending deformations. The periodic beat patterns are the result 
of a mechanical feedback that leads to self-organized bending waves along the axoneme. Using 
a theoretical framework to describe planar beating motion, we derive a nonlinear wave equation 
that describes the fundamental Fourier mode of the axonemal beat. We study the role of nonlin- 
earities and investigate how the amplitude of oscillations increases in the vicinity of an oscillatory 
instability. We furthermore present numerical solutions of the nonlinear wave equation for differ- 
ent boundary conditions. We find that the nonlinear waves are well approximated by the linearly 
unstable modes for amplitudes of beat patterns similar to those observed experimentally. 
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I. INTRODUCTION 



Cilia and flagella are hair-like appendages of eukaryotic cells exhibiting regular, wave-like 
oscillations [T]. Their ability to generate regular beat patterns plays an important role in 
many systems where motion on a cellular level is required [2]. Examples range from the 
propulsion of single cells, such as the swimming of sperm, to the transport of fluid along 
ciliated surfaces, such as the flow of mucus in the trachea. Ciliary and flagellar beat patterns 
are generated by an active structure called the axoneme which consists of nine microtubule 
doublets arranged in a cylindrical geometry [U El S]. A large number of dynein motor 
proteins are arranged between adjacent microtubule doublets and generate internal stresses 
within the axoneme that induce relative filament sliding and as a consequence axonemal 
bending [3 El H El El [ID] . 

Axonemal beat patterns have been the subject of several theoretical analyses trying to 
elucidate the mechanisms underlying the generation of regular beat patterns |TT1 [121 IISI 
[Ill[ISl[ISl[IZl[IHllISl[2nil2Il[221l2Sl[21]. Recent evidence suggests that the interplay of 
collectively operating motors together with the elastic microtubules constitutes a mechanical 
feedback that leads to oscillating instabilities [I3l [191 125]. The resulting traveling wave 
bending patterns can account for the experimentally observed beat patterns in bull sperm 

In the present article, we discuss the properties of self-organized beating patterns, ex- 
tending previous work in which the linearly unstable modes near an oscillatory instability 
were studied [191 • We present a nonlinear wave equation that describes the fundamental 
Fourier mode of planar axonemal beat patterns and derive analytically how the nonlin- 
earities determine the amplitude of the beat beyond the bifurcation point. Furthermore, 
we present numerical solutions of the nonlinear wave equation subject to three different 
boundary conditions. 

II. DYNAMIC EQUATIONS OF MOTION 

Motivated by the observation that the flagellar beat patterns of many sperm are approxi- 
mately planar we discuss the dynamics of the axoneme in a plane. Such planar beat patterns 
can be described by an effective, two-dimensional description of the three-dimensional ax- 
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onemal structure, in which the axoneme is represented by two elastic rods separated by a 
fixed distance a, corresponding to the axonemal diameter of 185 nm [121 123] ■ These rods are 
hnked by elastic structural elements and by active force generators, corresponding to the 
dynein motor proteins. To describe the relative motion of the two rods, we introduce the 
local sliding displacement A and the local shear force density / exerted by passive elastic 
and active elements, as illustrated in figure flj We denote by r(s,t) the two-dimensional 




FIG. 1: Schematic representation of the effective two-dimensional mechanics of planar beats with 
two elastic rods shding relative to each other due to the shear forces generated by active elements. 
Illustrated are the tangential shear forces /(s) and the local sliding displacement A(s). Elastic 
structural elements are indicated as springs. 

space curve parametrized by its arc length s describing the shape of the centre line of the 
axoneme of length L, at time t. As illustrated in figure |2| this shape can be characterized 
by the local tangent angle ip{s,t) such that 

r(s) = r(0)+ [ (cosV^(s'),sinV^(s'))ds' , (1) 
Jo 

where we have dropped the time dependence for notational convenience. In this two- 
dimensional geometry, the local sliding displacement and the local tangent angle are then 
related by 

A{s)=Ao + a{^{s)-tP{0)) , (2) 

where Aq denotes the relative sliding displacement at the base [101 EBl ES] . For simplicity, 
we ignore hydrodynamic interactions and describe the local hydrodynamic friction by in- 
troducing drag coefficients per unit length ^|| and C,± for movements in directions parallel 
and perpendicular to the axonemal axis, respectively. The dynamics of the axoneme with 
a bending rigidity k and an internal shear force density f{s,t) is then described by the 
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FIG. 2: Geometry of the flagellar deformation in the x,y-plane. The shape at a given time is 
described by the local tangent angle ip{s) as a function of the arc length s along the flagellum. 

following set of coupled nonlinear equations [19] 



r" - fi^'fr = a (/>' + /^") - n ((^")' + ^V") + ^ - , (4) 



where the primes denote derivatives with respect to the arc length s, i.e. ip' = dgip- The 
lateral tension r(s, t) ensures that the filament satisfies the local inextensibility constraint 
(r')^ = 1. Note that these equations can be derived from a full three-dimensional dynamic 
description of the axonemal cylinder, restricted to deformations in a plane and do not require 
the introduction of the effective, two-dimensional axoneme shown in figure [l] |2ll EZ] ■ 

III. BOUNDARY CONDITIONS 

The dynamic equations ([Sj-Q are complemented by boundary conditions. While the 
distal end (s = L) is typically free to move without external constraints, the basal end (s = 0) 
is subjected to external forces and torques. Matching the internal and external torques and 
forces at the ends determines the boundary conditions for the tangent angle ip{s,t) and the 
tension r(s, t) Different experimental conditions imply different boundary conditions at 
the basal end s = 0. Motivated by experiments in which the centre of the sperm head is held 
at a fixed position but potentially free to pivot [23], we describe the dynamics of the head 
angle '?/'(0, t) by introducing an angular elastic modulus kp and an angular friction coefficient 
7p. The set of general boundary conditions is summarized in table [Tj In this article we will 
discuss the cases of (i) a clamped head corresponding to the limit of large kp and (ii) a freely 
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pivoting head corresponding to fcp = and 7p = 0. 



At s = 


At s = L 


Kijj' + a f{s)ds — kpip 


- jpdti; = V' = 


Kip'" -af -ip'T = 


K^p" - a/ = 


kiIj'iI)" - afil;' + r' = 


r = 



TABLE I: Boundary conditions of sperm with fixed head position and free tail. We consider 
specifically the clamped head corresponding to the limit of large kp and the freely pivoting head 
limit with kp = and 7p = 0. 

In order to complete the description of the basal dynamics, it is necessary to specify the 
mechanical properties of the basal connection which determine the relative sliding between 
microtubules at the base flUl [26] . Recently, it has been shown that such basal sliding can 
have an important effect on the shape of the flagellar beat |23j . Following this previous work 
we characterize the visco-elastic coupling between microtubule doublets at the basal end by 
a basal elasticity ks and a basal friction 7^. The basal sliding displacement Ao(t) then obeys 
[23] 

jsdtAo = -fc.Ao - / fis)ds , (5) 
Jo 

and in the limit for large ks and 7^ basal sliding is suppressed. 



IV. OSCILLATORY DYNAMICS 
A. Fourier representation 

Time periodic beat patterns can be represented by the temporal Fourier modes ipn{s) of 
the tangent angle 

00 

^{S,t)= J2 ^nisy""^' . (6) 
?1=— 00 

The Fourier modes /„(s),A„(s) and r„(s) of the local shear force density f{s,t), the local 
shding displacement A(s,t) and the tension T{s,t) are defined identically. 

The motor proteins in the axoneme generate time dependent shear forces which induce 
dynamic sliding displacements A{s,t). The relation between sliding speed and force is a 
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collective property of the motors together with passive elements cross-linking the axoneme. 
This effective mechanical property of active and passive elements can be represented as a 
nonlinear relation in terms of the temporal Fourier modes [T^ 125] 

/i = «Ai + /3Ai|Ai|2 + 0(A5) . (7) 

The emergence of spontaneous oscillations is related to negative signs of the real and imag- 
inary parts of the linear response function a, resulting from the collective properties of 
many molecular motors coupled to an elastic element [191 ESI EH EE ED]. The collective 
effects arise from the dependence of motor transition rates on the state of the system, as for 
example, introduced by a load dependence of the motor detachment rate [HI [23] . 

B. Nonlinear waves 

The self-organized nonlinear dynamics of the axoneme can then be expressed by coupled 
differential equations for the discrete Fourier modes of the tangent angle and the tension. 
The experimentally observed beat patterns of sperm are dominated by their fundamental 
temporal Fourier mode, with higher harmonics contributing to less than 5% of the wave 
pattern [23]. In the following, we thus neglect higher harmonics of ip{t). 

To simplify the notation in the following, we drop the tilde when referring to temporal 
Fourier amplitudes, defining ip^s) = ■ipi{s),To = fo(s),r2 = ^2(5). We also introduce di- 
mensionless parameters u,a,(3 and Aq as defined in [a] Taking into account nonlinearities 
self-consistently up to cubic terms, equations ([3]),(|4]) and ([t]) lead to the following set of 
coupled nonlinear equations for the dominant modes 1/^, tq and f2, as rescaled in [X] 

d,{\ij\^) - 2a{^ + Ao - - «*(^* + AS - rmi^' + m-^o + rd,f2 

8% = 2-Re{ad,[{iP + Aq - ^(0))^*]} - dl{\i,\') + 2^ (l^l^Re {a} - Re{r^}) 

dlf2 = ad.liiP + Ao - ^(0))^'] - dS^) + |L [ai,^ - • (8) 

Here the dots denote derivatives with respect to the rescaled arc length s = s/L and complex 
conjugates are denoted by asterisks. In equation ([s]) we have also introduced the dimension- 
less linear contribution to the fundamental Fourier mode of the basal sliding displacement 
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Ao = ^ (m - [\is)d-s 

iu% + ks + a \ Jq 

The corresponding boundary conditions complementing equation ^ are summarized in |A] 

V. WAVE AMPLITUDES 

The above system exhibits an oscillatory instability or Hopf bifurcation, at which the 
unstable modes are described by a linearized wave equation [191 123] . In the oscillatory regime 
close to the bifurcation, finite amplitude solutions of the full nonlinear wave equation (|8| are 
similar to the linearly unstable modes. The growth of the amplitude and the changes of the 
shape of the beating mode with increasing distance from the bifurcation are determined by 
the nonlinear terms of equation (jsj). In the following we study the effects of nonlinearities 
near the bifurcation using a systematic expansion. 

Linearizing the nonlinear wave equation ([s]) in the limit of small amplitudes, the linearly 
unstable modes denoted by Uo{s) satisfy [121 [23] 

CMs) = , (9) 

subject to appropriate boundary conditions [IHl 123], where we have defined the linear oper- 
ator 

u) = iu; + - adl . (10) 

Note that the amplitude of the linear mode uq is not determined by the linear equation 
(9]). For convenience, we normalize uq such that |uo(s)|ds = 1. Equation (9]), together 
with the appropriate boundary conditions [TH |23], constitutes a boundary value problem. 
Nontrivial solutions exist only for pairs of critical values of the dimensionless frequency and 
response coefficient {(Xc^oJc) [IS1I2ZI- In the following we will denote Cc = C{ac,^c), where 
ac and Uc are the values of a and u at the bifurcation point. In figure [3} the line of critical 
values is indicated by the solid line representing C^Q clS db function of Uc- Note that there 
exists a discrete spectrum of such critical lines [T9] . 

We can express solutions of the full nonlinear problem for parameters a = ac + Sa and 
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uj = uJc + Su in the vicinity of the bifurcation point by an expansion of the form 



(11) 



Here e is a small dimensionless number that characterizes the distance from the bifurcation 
point by 

5a = pe'^e^ and 6iu = /ie^ , (12) 
where we have introduced the real coefficients p and fi as well as the phase 6. 






FIG. 3: Schematic diagram of the complex plane of motor impedance Oc- For a given dimensionless 
frequency iVc, there exists a critical value Oc describing an oscillatory instability. The line ac 
parametrized by cJc is shown in black. In the dashed region to the right of the line the system is 
quiescent, whereas to the left of (Xc it oscillates. The phase theta describes the orientation of a 
displacement Sa away from a bifurcation point in the complex plane such that 6a = \6a\e^^. 



Inserting the ansatz (11 ) into the wave equation nSv, we can solve this equation systemat 



ically order by order near a given bifurcation point. First, the linearly unstable mode Uq{s) 
is determined. Then with Uo{s) known, we can determine the static and dynamic tension 
profiles v(s) and w{s), in terms of uo{s) as detailed in[B| 

Matching terms to third order in e then leads to an equation for the nonlinear correction 
Ui{s) to the waveform 

CcUi = pe^^dluQ — ipuo — M{uo) , (13) 
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where the nonhnear terms A^('Uo) are given by 

^^{uo) =M + - mo(0))|mo + - mo(0)|'1 + Os {vu^ + wul) 



+ 



6 



5s(|%|'mo) - 2a(Mo + - Mo(0))|no|' 



(14) 



and Aq'^'' is the basal shding term evaluated at the bifurcation point, as defined by equation 



(B3) in the Appendix. 



Using equation ( 13 ) we can obtain a relation between the coefficients and 9 without 



calculating the nonlinear correction ui. This is achieved by introducing a function Uq ad- 
junct to the linear mode Uq which has the property C-cUq = and obeys CcUi = Z. 
The constant Z depends on 6,p,fi and is derived explicitly in O for three different mechan- 



ical conditions imposed at the basal end. Multiplication of (13) with Uq and subsequent 
integration leads to 

-pe'^ /" u^d'^uods + ifi f u^Uods + [ u^N'{\uo\^uo)ds + Z = . (15) 
Jo Jo Jo 

We can now discuss the emergence of the unstable mode ipls) and its frequency a) when 

starting at a bifurcation point at a = as illustrated in figure [3] moving from ac in a 



direction given by an angle 9, equation (15) describes the beating mode in the oscillatory 
region of the state diagram characterized hj a = ac + \Sa\e^^ . For a chosen value of 6, the 



values of p and p can be uniquely determined from the complex equation (15). Equation 



(12) then describes the increase of the amplitude as 



P 

while the frequency changes at the same time by 

60 = — \6a\ 
P 



(16) 



(17) 



Note that these behaviors depend on the angle 6 chosen. Two special situations are 
of interest. There exists in general a specific choice 6 for which /i = 0, i.e. in this case 
the frequency does not change when moving in the corresponding direction away from the 
bifurcation line. Examples for such lines in the complex plane along which the frequency of 
the unstable modes remains the same as at the bifurcation point dc are displayed in figure 



[4| A second special choice 9 = 9\\ is the direction tangential to the bifurcation line q;c(c<Jc)- 



For this choice, equation (15) becomes singular with p ^ oo and p/yU = Idac/dcud. In this 
case the amplitude e remains zero, but the frequency changes along the bifurcation line [T^ . 
Note also, that the shape and frequency of the beating mode at a point a do not depend on 
the reference bifurcation point from which it is reached. 

The above method permits us to calculate the amplitude and frequency of beating modes 
close to the bifurcation line. In order to study the influence of nonlinearities on the shapes 
of these modes, we make use of the above analytical result to solve the nonlinear equations 
numerically. 

VI. NUMERICAL SOLUTIONS TO THE NONLINEAR WAVE EQUATIONS 

Periodic and planar beating patterns are solutions to the nonlinear wave equation (|8 



together with the boundary conditions given in table III which constitute a boundary value 
problem, that can be solved numerically by a shooting and matching procedure [31]. Note 
that the wave equation is invariant with respect to changes in the overall phase of ip{s). To 
remove this degeneracy we impose an arbitrary condition for the phase of iIj{s) at the base 
s = 0. 

In order to obtain numerical solutions that satisfy the wave equation with given boundary 
conditions, we first determine an approximate solution tp{s) ~ euo(s),fo(s) ~ e^f (s),f2(s) ~ 
e^w(s) close to a bifurcation point cuc, making use of the method discussed in the preceding 
section. This allows us to determine the amplitude e and the functions v{s),w{s) for a 
solution close to the bifurcation point ctc^c- This approximate solution is used as seed for 
the shooting and matching procedure to solve equation Q. 

With this procedure, we generate a sequence of numerical solutions starting from the 
bifurcation line and moving in a direction in which the frequency ui remains constant. Ex- 
amples of these solutions to equation (|8| are displayed in figure |4] for different boundary 
conditions and starting from different bifurcation points. 

Parameter values used in these calculations are = 3.4 x 10~'^N ■ s ■ m~^, = C,±/2, 
K = 1.7 X 10~^^ N ■ m^, L = 58.3 /im as estimated for bull sperm fiagella E2]- The value 
of f3 has not been measured. We choose /3 = 42 for which the bifurcation is supercritical. 

The region of stability of the non-oscillating state is indicated in figure |4]^ A, C,E). The 



10 



oscillatory instability occurs along the solid black line. The real and imaginary parts of 
for these solutions are displayed in figure |4]^B,D,F) as a function of the dimensionless arc 
length s for distinct values of a, located along the grey line shown in figure |4]^A,C,E). As 
indicated by the insets, the amplitude A = {^(s) \ds of the modes grows continuously with 
increasing distance from the instability. In the limit of small amplitudes, it obeys equation 



(jiej) with e ~ A. 

In contrast to the amplitude, the shape of the beat patterns changes only weakly as illus- 
trated by figure [5] in the Appendix. Solutions to the linearized equations therefore provide 
good approximations to the full nonlinear problem in the range of examined parameters. 



VII. CONCLUSIONS AND OUTLOOK 



The beating of cilia and flagella in viscous media is the result of the co-operative action 
of many dynein motors interacting with the elastic microtubules within the axoneme. The 
solutions to the linearized dynamic equations [19] have been shown to be good approxima- 
tions for experimentally observed planar beat patterns of bull sperm [23j. However, the 
flagellar beat is a fundamentally nonlinear problem and cannot be satisfactorily described 
by a linear theory. 

Here, we have presented a theory of flagellar beat patterns taking into account the leading 
nonlinearity. We derived the nonlinear wave equations for the fundamental modes of the 
beat shape and the lateral tension ([s]). Furthermore, by solving the nonlinear wave equation 
numerically we explored how the shape of beat patterns changes as the beat amplitude 
increases. We found that for the parameter values and amplitude range examined here, 
the shape of the beat patterns of flnite amplitude remain qualitatively similar to the linear 
waveforms obtained near the instability in the limit of inflnitesimally small amplitudes. 
Experimentally observed wave forms of beating sperm can therefore be approximated by 
linear modes of ad hoc amplitude [23]. Existing nonlinear discussions of slender rod dynamics 
in the viscous regime have focussed on the dynamics of passive rods [331 EH |35] , or described 
actively beating fllaments subject to pre-described internal force distributions [SH]- Our 
analytical results provide the framework to discuss the lowest order nonlinear effects of self- 
organized axonemal beat patterns allowing us to bridge the gap between the linear and 
nonlinear regime without having to resort to a full simulation of the system. 
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FIG. 4: Examples of beating modes, which are solutions to the wave equation (|8| for different 
boundary conditions. The location of these solutions is indicated in the complex plane of values 
of the linear response function of motors a. The chosen examples correspond to a situation where 
starting from a Hopf bifurcation, the amplitude of the unstable mode increases when moving along 
a path in the a plane for which the oscillation frequency is constant. (A) Solutions for clamped 
head with basal sliding with frequency {jjjlTx ~ 26 Hz in the a plane (grey line). The parameters 
kg and 7^ were chosen as determined in [23^ such that the beat patterns shown in (B) resemble 
experimentally observed ones. The amplitude j4 as a function of |5a| is displayed in the inset. 
(B) Real and imaginary part of the waveform ^/^ as a function of dimensionless arc length s for 
different points along the grey line in (A). (C,D) Same as (A,B) but for the clamped head boundary 
conditions without basal sliding and frequency ~ 28 Hz. Here we choose the bifurcation on the 
first branch of unstable modes. (E,F) Same as (C,D) but for a freely pivoting head without basal 
sliding and frequency ~ 5 Hz. 
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Our results are independent of the specific molecular mechanisms underlying the collective 
action of motors. However, our theory does therefore not allow us to predict the influence 
of experimentally controlled parameters such as ATP concentration or temperature on the 
beat shape. To understand the influence of such parameters on the beat, more specific 
models of motor action are required. The oscillation frequency of the system is selected via 
the self-organization of dynein motors and microtubules. This frequency selection involves 
the full frequency dependent impedance of motors and thus also depends on details of 
the underlying mechanisms [12]. A key open challenge for the understanding of flagellar 
dynamics is therefore to understand the parameters which govern frequency selection based 
on the properties of dynein motors in the system. Also, there exist several branches in 
parameter space along which linear instabilities occur. Which of these branches becomes 
unstable first and governs the behavior of the flagellar beat also depends on molecular details 
and will be an important challenge for future work. 
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APPENDIX A: DIMENSIONLESS PARAMETERS AND BOUNDARY CONDI- 
TIONS 



Aq = Ao/a OJ 
fi{s) = ^Ti{s) ks 

TABLE II: Summary of the relations between the physical parameters and the dimensionless quan- 
tities. 

Table [H] provides the definitions of dimensionless parameters used in this work. The 
mechanical conditions at the boundary (see table B impose boundary conditions for the 
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fundamental modes ipi^)} ''"o(s) and T2(s) as summarized in table III, where we have intro- 
duced the nonlinear contribution to the dimensionless basal sliding displacement 



A, 



(nl) 



Ao 



iuj^ + ks + a Jo 



7/^(0) + Ao|'(V^(s)-^(0) + Ao)ds . 



(fep + iw7p)^(0) = ^-(0) + a [ [V'(s) + Ao - ^(0)] ds 

Jo 

(s) + Ao - ^(0)|2(^(s) + Ao - V(0))] ds 



Jo 

= a^(0) + V'(0)to(0) + ri0)f2{0) + ^(2|AopV'(0) + (Aq) V(0)) 

dfo 

ds 
df2 



5=0 



-55(|V'(0)|2)+2Re{QAor(0)} 



-V'(O)'0(O) + aAoV'(O) 



5=0 

^(1) = 



v^(i) = am) + Aj'^') + v^(o)) + pmi) + Ao - m\''m) + ao - m) 

T0(1) = 
T2(l) = 



TABLE III: The boundary conditions for the fundamental modes il^{s), to(s) and T2(s). 

Note that for the limiting cases under consideration in the main text, the boundary 
conditions simplify as follows. In the general clamped head case we have ^^(O) = 0, in the 
absence of basal sliding we furthermore have Ao = 0, and Aq"''' = 0. For a freely pivoting 
head without basal sliding, kp = 0, 7p = 0, Ao = 0, and Aq^^^ = 0. 

APPENDIX B: NONLINEAR PERTURBATION CALCULATION TO SECOND 
ORDER 

Formally expanding the linear operator C close to the bifurcation point, we write 

ij) = C,- 5adl + i5uj . (Bl) 



Substituting the ansatz of (11 ) and (12) into the nonlinear wave equation m\ reproduces to 



linear order the equation describing the linearly unstable modes (|9| supplemented by the 
appropriate boundary conditions (explicitly described in |27j ) . 
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Matching terms to second order in e then leads to 

v{s) = 2Re{a,d-s[{uo + Af^ - no(0))M*]} - dl{\uo\^) + 2^ {\uo\^Re {a,} - Re{n*Mo}) 

c 

iu{s) = acdsliuo + ^ - mo(0))uo] - ds{uoUo) + (ac(^io)^ - uqUo) (B2) 
where 



(c) ttn 



AJ,^^ = -7^— - / Ms)drs , (B3) 

IWcTs + Ks + Jo 

is the amphtude of basal sliding to linear order. The above system of equations, together with 
the appropriate boundary conditions p7], allows us to obtain v{s),w{s) for given solutions 
Uq{s) of the linear problem. 



APPENDIX C: NONLINEAR PERTURBATION CALCULATION TO THIRD 
ORDER 

We define a function Mq adjunct to Uq such that the integral 

"1 /•! 



Jo Jo 



1 

ui{iuJcU^ + "Uq - acUQ 







is independent of ui. In order to eliminate ui{s) from the bulk term in the above expression 
we require that CcUq{s) = 0, implying that Uq{s) satisfies the same differential equation as 
the linear modes Uo{s). The condition that terms depending on ui vanish at the boundaries 
s = and s = 1, then defines the boundary conditions for Uq. The so-defined function 
Uq then leads to J^UqCcUi = Z, where the complex constant Z = Z{p,fx,9) does not 
depend on the unknown correction function ui. This complex constant characterizes the 
behavior of the system away from the bifurcation in the direction of 6 and can be determined 
from the linearly unstable mode Uq and its adjunct function Uq only. In the following we 
state explicitly the boundary conditions specifying the adjunct function and determine the 
constant Z for the various mechanical conditions of interest (clamped head without basal 
sliding, clamped head with basal sliding, freely pivoting head without basal sliding). These 



results were obtained by matching terms of 0{e)-0{e^) of the perturbation ansatz (11) 
substituted into the nonlinear wave equation 
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1. Clamped head without basal sliding 

The boundary conditions for the conjugate hnear solutions W^ls) are given by 
M+(0)=0, <(0) = 0, M+(1) = 0, Mo+(l) = , 

which leads to 



/ u+C,ui = -u+{0)A-u+{l)B = Z 
Jo 



where we have introduced 



5 = pe%(l) + ^|Mo(l)|'Mo(l) . 



2. Clamped head with basal sliding 



(CI) 



The boundary conditions for the conjugate linear solutions Uq{s) are given by 

1 

7i+(0)=0, <(0)=0, n+(l) = 7^7^— ^— -4(1), ^■o^(l) = 0, 



(C2) 



which leads to Z of the same form as (C2), but with a more complicated expression for A 
and B which are now given by 

A =pe%(0) +^(2|Af)|2Mo(0) + (A^^))^^^^) + t;(0)Mo(0) + W7(0)<(0) 



(c) 



-etc. 



+ 



_ , , ^ . pe%(0)-A + / Ar(|wo(^)pwo(s))ds 
iWcTs + ks + ac luJc V Jo 

1 



+ 7 



iuc% + ks + ac 



Af^(pe'^(icje7s + ks) - ifi%ac) - P^c I \uo{s) + A^'f (mo(s) + A'^^>)ds 



(c)|2/ 



,(c)^ 



with A/" as defined in ( 14 ) 



3. Freely pivoting head without basal sliding 



The boundary conditions for the conjugate linear solutions are given by 



M+(0) = 0, u + il) = 0, <(1) + i-4(0) = 



UJc. 



u + {0) - ae«(0) + n+(0) - <(1)) = 



(C3) 
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which leads to 



UqCcUi -- 
where we have introduced 



LUr 



A =pe%(0) + t;(0)Mo(0) + u'(0)mS(0) 

B =^iio(l) + |wo(1)|'mo(1) 

otc \ac\ otc 



pe 



id 



C =^M0) - P / Ms) - uo{0)\\uois) - uo(0)) 

ttc Jo 



ds 



etc 



+ - 



1 



pe%(0) + i/i(Mo(0) - — ito(O)) + / ^i\uo\^Uo)ds 

ttc Jo 



(C4) 



APPENDIX D: SHAPE CHANGES OF THE NONLINEAR SOLUTIONS 
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FIG. 5: Examples of the shape of the nonlinear solutions as they grow from the linearly unstable 
modes, corresponding to the 26 Hz solution of the second branch of unstable modes in the case of 
clamped head boundary conditions with basal sliding. Shown are the real and imaginary parts of 
the rescaled nonlinear solutions with the lighter curves corresponding to larger amplitudes. The 
amplitudes A (in rad) of the selected modes are indicated by black dots in the diagram that shows 
how the amplitude grows as a function of the distance \6a\ to the bifurcation. 
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